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The renormalization group method is a successive integration over the fluc- 

CN ' tuations which are ordered according to their length scale, a parameter in the 

Qv ' external space. A different procedure is described, where the fluctuations are 

0^ . treated in a successive manner, as well, but their order is given by an internal 

^^ , space parameter, their amplitude. The differential version of the renormal- 

I , ization group equation is given which is the functional generalization of the 

Q I Callan-Symanzik equation in one special case and resums the loop expansion 

i-G ' in another one. 
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I. INTRODUCTION 

The renormalization group method is usually applied in two different manners. The 
original way is to provide an insight into the scale dependence of the coupling constants [|l|. 
Another, more recent use is to perform a partial resummation of the perturbation expansion 
by making an infinitesimal change of the cutoff in a time and using the functional formalism 
0" 101 • Both goals are realized by the blocking procedure, the successive elimination of the 
degrees of freedom which lie above the running ultraviolet cutoff. The resulting evolution 
equation yields the dependence of the coupling constants in the cutoff which is a scale 
parameter introduced in the space-time, the external space. A difficult problem in this 
program is the lack of explict gauge invariance. Since the separation of the modes according 
to their length scales is not gauge invariant, the gauge symmetry is unavoidably lost during 
the blocking. 

We propose here a different slicing procedure of the path integral, the gradual increase of 
the amplitude of the fluctuations. We introduce a cutoff in the internal spaceQ by constrain- 
ing the amplitude of the fluctuating field and obtain the evolution equation by the gradual 
release of this constraint. Since the field variable is usually dimensional this procedure can 
be interpreted in terms of a scale dependence, as well. More precisely, the mass parameter 
of the lagrangian serves as a cutoff for both spaces thereby providing a dynamical corre- 
spondence between the scales in the internal and the external spaces. An important feature 
of this scheme is that the scale parameter, the mass in particular, is introduced by hand 
in the internal space. The length scale of the external space is, on the contrary, induced 
dynamically. Such a dynamically generated length scale in the external space appears more 
natural and avoids the problems mentioned above. In fact, having imposed the cutoff in 
the internal space the gradient expansion, being an approximation procedure in the external 
space, goes through without problem. The gauge invariance can be maintained, as well, the 
only requirement is that the suppression of the fluctuations must be achieved in a gauge 
independent fashion, a condition easy to comply with in the internal space. This latter 
improvement is somehow reminiscent of the achievement of the Higgs mechanism where a 
dynamically generated mass is produced without spoiling the gauge symmetry. 

There are further advantages in making in the internal space the analogue procedure of 
the blocking. One is to get a new insight into the origin of the anomalous dimension. To see 
this we recall that the internal space has actually been used for the purpose of establishing 
the renormalization group. The renormalization conditions was imposed at a non-vanishing 
background field field in ref. 0] providing the scale for the internal space. The resulting 
renormalization group function was computed in the one-loop level and found to be identical 
with the usual functions, resulting from the renormalization conditions introduced in the 
external space. This agreement between the scale dependence in the external and internal 
space is not generally present and the two-loops results should already yield a difference. 
In fact , the scales inferred from these two spaces are related by the terms of the lagrangian 
which mix the internal and the external spaces, the terms with space-time derivatives. The 



^The internal and the external spaces stand for the space of the field amplitude and the space-time 
manifold, respectively. 



simplest of them is the wave function renormahzation constant, (J)b{,x) = Z^/^</)/j(x), which 
is trivial for a scalar model in the one-loop approximation. The anomalous dimension is 
nothing but the demonstration of the difference between the scale dependence in the internal 
and the external spaces. 

Another bonus of the internal space renormahzation group is the generalization of the 
Callan-Symanzik equation for the functional formalism. The Callan-Symanzik equation de- 
scribes the dependence of the Green functions on the mass 0. Since the mass parameter 
controls the width of the peak in the wave functional around the vacuum configuration the 
mass can be thought as the simplest cutoff parameter in the internal space. The similarity 
of the renormahzation group coefficient functions obtained by the bona fide renormahzation 
group and the Callan-Symanzik equation demonstrates that the length scale of the exter- 
nal space induced by the mass can be identified by the cutoff. The functional evolution 
equation resulting from the infinitesimal change of the mass is the generalization of the 
Callan-Symanzik equation without relying on the multiplicative renormahzation scheme. 
Such an extension, presented below can be used to study non-renormalizable models or 
effects close to the cutoff. 

When the procedure in the internal space is performed by the help of the bare action 
the resulting flow provides a resummation of the loop expansion. This scheme is specially 
advantegous for models with inhomogeneous saddle points because the blocking does not 
distorts the tree level structure. 

The renormahzation group step can be expressed by means of the bare action or the 
effective action. The two schemes become equivalent in the infrared limit. This can be 
seen easily by remarking that the loop contributions are suppressed in the infrared limit, 
where the ultraviolet cutoff tends to zero with the number of degrees of freedom. Since the 
difference between the bare and the effective action comes from the loop contributions the 
gradient expansion produces the same bare and effective action in the vicinity of the infrared 
flxed point. 

The different standard methods to obtain the renormahzation group are briefly discussed 
in Section II. The internal space evolution equation, our new renormahzation group scheme 
is introduced in Section III. Section IV contains the application of the general strategy of the 
internal space renormahzation to obtain the generalization of the Callan-Symanzik equation. 
The possibility of resumming the loop expansion by solving the evolution equation is shown 
in Section V. The Section VI is for the summary. The appendices give details about the 
Legendre transformation and the gradient expansion. 

II. RENORMALIZATION GROUP SCHEMES 

In this Section we briefly recapitulate some methods the renormahzation group equation 
can be obtained. 

The traditional fleld theoretical methods for the renormahzation group equation are 
based on the simpliflcation offered by placing the ultraviolet cutoff far away from the scale 
of the observables. Such a separation of the scales removes the non-universal pieces of the 
renormalized action and the rather complicated blocking step can be simplifled by retaining 
the renormalizable coupling constants only. The underlying formalism is the renormalized 
perturbation expansion, in particular the multiplicative renormahzation scheme. The usual 



perturbative proof of the renormalizability asserts that the renormahzed field and the Green 
functions can be written in terms of the bare quantities as 

Gnipi,---,Pn]gR{fJ'))R = 2^2 igR,gB,-j Gn{pi, ■ ■ ■ , Pn] 93, A) B ( ^ + *^ ( ^2 ) ) ' ^^^ 

where A is the cutoff and the renormahzed couphng constants are defined by some renor- 
malization conditions imposed at p"^ = ^^. The evolution equation for the bare and the 
renormahzed coupling constants result from the requirements 

±.0. . ±Z-,0. . 0. (2) 

±Gb = ^ZiG^ = 0. (3) 

djji djj, 

Note that the non-renormalizable operators can not be treated in this fashion because the 
0{p/A) contributions are neglected in (|l|). The renormalization of composite operators 
and the corresponding operator mixing requires the introduction of additional terms in the 
lagrangian. Another aspect of this shortcoming is that these methods are useful for the 
study of the ultraviolet scaling laws only. The study of the infrared scaling or models where 
there are several non-trivial scaling regimes |1^ require the more powerful functional form, 



introduced below. 

The third conventional procedure is the Callan-Symanzik equation which is based on the 
change of the bare mass parameter, 

where Z^2 is the renormalization constant for the composite operator 0^(x) and C^^^^p is 
the Green function with an additional insertion of 0^(p = 0). One can convert the mass 
dependence inferred from the Callan-Symanzik equation into the momentum dependence by 
means of dimensional analysis and the resulting expression is usually called a renormalization 
group equation. 

The functional generalizations of the renormalization group method which are based 
on the infinitesimally small change of the cutoff allows up to follow the mixing of non- 
renormalizable operators, as well, and to trace the evolution close to the cutoff. Another 
advantage of these methods is that the renormalization group equation is either exact or 
holds in every order of the loop expansion. 

The traditional blocking [|^] in momentum space yields the Wegner-Houghton equation 
which is based on sharp cutoff in the momentum space. Consider the action with cutoff 
k in the derivative expansion, 

oo „ 

S, = J2 d'xUi-\<P{x),d,] (5) 

n=0 

where L^^" [0(x), 9^] is an n-th order homogeneous polynom in the gradient operator d^ and 
a local functional of the field variable. We separate the quadratic part of the action in the 
field. 
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2 
where 

to organize the loop expansion. The blocking step consists of the elimination of the Fourier 
components of the field 4>{x) within the shell k — Ak < p < k. Let us write 

0(x) = 0(x) + 0'(x) (8) 

where the first and the second term in non- vanishing in the Fourier space ioi p < k — Ak 
and k — Ak < p < k, respectively. The blocked action can be written as 

^ ^-Skl4>+'4>]-htrG^^[4>]+- (9) 

where ip is the saddle point, the trace is within the functional space spanned by the plane 
waves k — Ak < p < k, and the dots stand for the higher loop contributions. The density 
of the modes to be eliminated, Q^^kk'^'^ , Q^ being the solid angle, is chosen to be small. 
This makes the volume of the loop integral appearing in the computation of the trace small, 
too. One expects that the n-th loop contributions which contain an n-fold integration in 
this small volume will be suppressed as Ak — ^ since the integrand is non-singular in the 
integration domain. A more careful argument shows that the small parameter is actually 
the ratio of the modes to be eliminated and left in the system 0, [|T0[ assuming that the 
amplitude of each Fourier mode is of the same order of magnitude^. With this condition 
kept in mind the higher loop contributions are dropped yielding the evolution equation 

S,[^] - 5.-A.M = -itrG-i[0] (l + O (^)) (10) 

for vanishing saddle point, ip = 0. All we needed in deriving this equation was the availability 
of the loop-expansion for the blocking step (P). The simplification, the suppression of the 
higher loop contribution, results from the appearance of a new small parameter Ak/k and 



(p!0|) is valid for weakly coupled models only. 

We may call this the scheme the functional form of the bare renormalization group, 
(0), since the resulting renormalized trajectory is in the space of the cutoff theories. The 
truncation of the gradient expansion 



^It is enough to introduce an infrared cutoff and consider the elimination of the modes one-by- 
one when this latter condition fails to become convinced that the limit Ak — > can not always 
suppresses the higher loop contributions to the blocking relation. Such a problem shows up when 



there is a non-trivial saddle point to the blocking which enhances certain modes |12]. 



Sk = I d'^x 



(d,^)\x) + U,i<Pix)) 



(11) 



simplifies the functional equation (1l0|) into a partial differential equation 






(k^ + dlUu{p))U + -dMp) 



N-l 



(12) 



for an iV-component field of modulus p. 

The inclusion of the 0{d^) higher order terms in the action is problematic because the 
sharp momentum space cutoff generates non-local interactions |11|. 

A similar method was developed for the effective action [^]- |^ where one introduces the 
cutoff parameter k in the quadratic part of the action G^^{p) — > G^^{p), with 

G-,\p) = f[l)G-\p), (13) 

where /(k) approaches 1 and oo for k = oo and k = 0, respectively. We have hope to 
generate local, gradient expandable interactions only if /(k) is an analytic function. Thus 
locality almost excludes the complete suppression of the modes, the choice /(/t) = oo for 
K, > kq < oo. The exception to this rule consists of the periodic, analytic functions, 

G^\p + K) = G-,\p) (14) 

where the components of the vector K are integer times k. An example is the lattice 
regularization where the complete suppression of the modes outside the first Brillouin-zone 
can be achieved without generating non-local interactions. When the modes are suppressed 
gradually one characterizes the suppression functions by the condition that it changes the 
most within the interval 1 — e<fi:<l-|-e. The limit e — *> is called the sharp cutoff. The 
path integration consists of the elimination of the modes p > k and the evolution equation 

dk /'z}[(/.]e"^'=[^l+i^'-'^- = - f D[(j)]dkSk[(j)]e-'^''^^^^^^''"^'' (15) 

is obtained for the generator functional. This relation is finally converted into an equation 
for the effective action. Since the effective action for the modes p > k is followed in this 
manner, these methods appear the analogues of (0), the evolution equation for the running 
coupling constants. 

The application of ( [T5| ) is a safe non-perturbative step only for a differentiable suppression 
function /(k). When the /(k) = oo for finite values of k or when /(k) has infinite derivative 
(say e — s> 0) then ([T5| ) is invalid. Instead one has to consider the finite difference equation 
k -^ k — Ak and the corresponding functional Taylor expansion brings back the assumption 
concerning the applicability of the loop-expansion. Thus (|T5|) produces an exact relation for 
smoothly suppressed modes only. 

The internal space renormalization group method presented in the next section gen- 
eralises the third traditional way of introducing the renormalization group, the Callan- 
Symanzyk equation by controlling the amplitude of the fluctuations in the internal space 
only. The coresponding suppression function can be choose as smooth as possible, /(k) = A, 
in order to suppress the higher order derivative terms generated by the blocking. 
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III. EVOLUTION EQUATION 

Our goal is to obtain the effective action r[0] of the Euchdean model defined by the 
action S'b[0], by reducing the renormalization group strategy into an algorithm to solve the 
theory. The connection to the external scale dependence will be considered later. The usual 
Legendre transformation (c.f. Appendix A.) yields 

gW^W^ /'^[0]g-5B[<^]+lj.</'. (16) 

and 

W[j] + T[<P]= 1 3A.=3-4>. (17) 

the source j is supposed to be expressed in terms of 

^ 5W\3] 



Sjx 



(18) 



A cutoff A is assumed implicitly in the path integral and S'b[0] stands for the bare, cutoff 
action. 

A sharp cutoff could be introduced in the internal space as 

Z^ = l[ /rf0,e-^[^l, (19) 

but it leads to unwanted complications and will be replaced by a smooth cutoff as follows. 
We modify the bare action 

Sb [0] -^ Sx [0] = \Ss [0] + Sb [0] (20) 

in such a manner that the model with A ^ oo be soluble. This is achieved by requiring 
that the minimum (p^o of the local functional Ss[(/)] > 0, is nondegenerate. In fact, the path 
integral for A = oo contains no fluctuations. The role of the new piece in the action is to 
suppress the fluctuations around 0oo for large value of A and render the model perturbative. 
The control of the amplitude of the fluctuations offered by the parameter A is considered 
as a smooth cutoff in the internal space. The difference between our strategy and the other 
methods is that the internal space cutoff appears multiplicatively in the action and influ- 
ences all modes simultaneously while the external space cutoff is modifying the momentum 
dependence. In other words, it is the order the different fluctuations are treated as the cutoff 
is decreased what distinguishes the internal and the external space renormalization group 
schemes. 

We plan to follow the A dependence of the effective action by integrating out the func- 
tional differential equation 

dxT = J^x[r] (21) 

from the initial condition 






knitSs[(j)] + Sb[ 



(22) 



imposed at Xinu ~ oo to A = 0. (^iD can be interpreted as a generalization of the Callan- 
Symantzik equation becaue both generate a one-parameter family of different theories orga- 
nized according to the strength of the quantum fluctuations 0. So long as the parameter A 
introduces a renormalization scale, /i(A), the trajectory rA(;i)[0] in the effective action space 
can be thought as a renormalized trajectory. Another way to interpret (^Tj) is to consider 
its integration as a method which builds up the fluctuations of the model with A = by 
summing up the effects of increasing the fluctuation strength infinitesimally, A ^ A — AA. 
Notice that the gauge invariance of the evolution equation ( PT| ) is obvious when Sg is gaue 
invariant. The gradient expansion is compatible with (|2T| ) if the suppression is sufficiently 
smooth in the momentum spcae, i.e. 5*^ is a local functional. 
The starting point to find J-'\[T] is the relation 



d,r\ 



-dxW[j] 



5j 



■dxj + dxj.^ 



-dxW[j], 



(23) 



A and being the independent variables. This relation will be used together with 



dxW[j] 



-W[j] 



D[(P]Ss[(P]e 



-XSs 



-W[j]g 
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oW[j] 



It is useful to perform the replacement 

r[0] -^ \s. 

which results the evolution equation 

dxT[<p]=e-'^^^^S^ 

The next question is the choice of Ss [(f> 
term, 



S_ 



+ r[0] 



,W[j] 



S.. 



(24) 



(25) 



(26) 



5. 



^xM^ 



The simplest is to use a quadratic suppression 

(27) 



x,y 



yry 



1 

2* 



M- 



If a gauge symmetry should be preserved then such an 5*^ 

SM = Sb[<I>] 



is not acceptable and 



(28) 



is the most natural choice. The evolution equation (|2q) then sums up the loop expansion 
and produces the dependence in h. We return now to the case of a simple scalar field without 
local symmetry, (|27|) . The corresponding evolution equation can be ontained from (|26|) , and 
considering the relation (|A4|) between the functional derivatives of W[j] and r[0]. 



^Note that the inverse mass is proportional to the amplitude of the fluctuations. 



^ Jx,y '- 



rM 



x,y 



x,yry 



M 



x,y 



rS + >^^^,y 



2 Jx,y 

where the functional derivatives are denoted by 



-1 



Yin) 



(29) 



(30) 



^Xl 



(^) reads in an operator notation 



-TilM 



XM + r(2) 



-1 



(31) 



We should bear in mind that T}^' ^ remains a functional of the field 



It is illuminating to compare this result with the evolution equations presented in refs. 
• what one can obtain by means of (p!3D and (p6|), 



d,T\ 



2Tr \ d,G-,' ■ 



Gz' + r(2) 



(32) 



The role of Ss is played here by the propagator Gk{p) which contains the external space 
scale parameter k to control the supression of the fluctuations. The formal similarity with 
(PSP reflects that the different schemes agree in "turning on" the fluctuations in infinitesimal 
steps. But the internal space scheme operates with a suppression term which is regular in 
or even independent of the momentum. 

The evolution equation can be converted into a more treatable form by the means of the 
gradient expansion, 



where the notation /^ = /(0^.) was introduced. This ansatz gives 



(33) 



(34) 



f(2) 






where the /*■"■' (0) = (9^/(0). Such an expansion is unsuitable for iy[j] due to the strong 
non-locality of the propagator but might be more successful for the effective action where 
the one-particle irreducible structure and the removal of the propagator at the external legs 
of the contributing diagrams strongly reduce the non-local effects. The replacement of this 
ansatz into ( p6D gives (c.f. Appendix B.) 



dxUM = lJ 



Mip) 



p\M{p) + Zxi^)p^ + U?{(P) 



d,z,{<p)=y^M{p) 

+ 2Z«(0) 



.(2) 



^r(</>) 



AA<(p) + Z,(0)p2 + ^7f (0) 
(, ,^2(z«(0y + f/f (0)) +Z«(0y/rf 



(Z«(0y + f/f (0))' (AnA^(p) + 2Z,(0)) 



-d^ 



(1) 



{XM{p) + Zx{(P)p^ 
(0)(zW(0y + [/f(0) 



f/f) 



(Ap^g^A<(p) + 2ZA(0y) 
AA<(p) + ZA(0)p2 + f/f (0) 



+ ^ 



2 (Zf (0)/ + f/f (0))' (A9^A^(p) + 2Z,(0)p, 



A^(p) + Z,(0)p2 + [/f(0) 



(35) 



where Jp = J ij^ and we assumed that d^M.{p) is proportional to p^. Since the order the 
fluctuations are treated is different for the internal and the external space methods, (|35D is 
different, as well, from the evolution equations of |0]- 0]. The final solution, correpsonding 
to A = in our case and fc = in the others, is the effective potential expressed in terms of 
the bare and the renormalized coupling constants, respectively. 



IV. MASS DEPENDENCE 



The simplest choice is A = m^ with 



•'^^x.ii ^x^ 



(36) 



which minimizes strength of the higher order derivative terms generated during the evolution 
by being a momentum independent suppression mechanism. The evolution equation is the 
functional differential renormalization group version of the Callan-Symanzik equation. 



d^2V\ 



iTr 
2 



m 



5 + r(^) 

'Jx.V \ ^ Til 



'x,y 



x,y 



(37) 



The projection of this functional equation onto the gradient expansion ansatz gives 

1 /• 1 



a„2f/(0) = - 

a^.Z(0) = i| 



2ypZ(0)p2 + m2 + f/(2)(0) 

Z(2)(0) 



(Z(0)p2 + m2 + [/(2)(0))2 

,^, p'/dZ^^\<P) + 2 (Z(i)(0)p2 + ?7(3)(0) 
2Z«(0)— — ^ —^ ^— - 

(Z(0)p2 + m2 + f/(2)(0))^ 



2Z(0) 



Z«(0)p2 + f/(3)(0) 

(Z(0)p2 + m2 + f/(2)(0))' 
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■ip 



_^Z(0)Z«(0) 



d 



+ ¥^^(0) 



(Z«(0y + f/(3) 



ci 



Z«(0y + f/(3)(0) 



(38) 



It is important to bear in mind that we are dealing here with a well regulated theory and that 
the procedure described here does not aim at removing the external space cutoff A which 
remains an important external parameter. In the usual, external space renormalization 
schemes the blocking provides the cutoff and the regulator for the models. In our case the 
cutoff in the internal space can not replace the usual external space cutoff A and this latter 
must properly be implemented from the very beginning, c.f. the remark after (|18|). In this 
case the loop integrals occuring in the evolution equations are finite. 

Let us now simplify the differential equation for f/(0) and Z{(j)) by integrating over p in 
(^) with sharp momentum cutoff A0 in four dimensions. 



d^2Z{(j)) 



1 


327r2Z(0) 
1 


327r2Z(0) 



A^ 



m2 + [/(2)(0) 

m 



In 1 + 



Z(0)A2 



m2 + t/(2)(0)_ 



Z2(0) V2 



ZW(0))'-Z(0)Z(^)(0))ln(l + 



Z(0)A2 



m2 + [/(2)(0)_ 



+ 



43 



zi<p)z^'\<p)-^{z('\<p) 



Z2(0)V~'"'~ '"' 12 

1 Z«(0)f/(3)(0) 1 (f/(3)(0) 



Z{<P) (m2 + f/(2) (0)) 6 (m2 + f/(2) (0))^ 



(39) 



In the approximation Z = 1 we obtain 



dm^U{(f)) 



m' 



+ U^^H^) 



327r2 



In 1 + 



A^ 



m2 + [/{2)(0)_ 



(40) 



after removing a field independent term. In order to simplify the scaling relations we consider 
the regime rri^ ^ f/^^-*, where 



a„2 [/((/)) 



327r' 



1 ,„ ('!!!!+4! V^'W 



m^ 



(41) 



The asymptotic scaling formula of the bare renormalization group (0) for d = Ain the same 
regime is 



^ The sharp momentum space cutoff generates nonlocal interactions. Since these nonlocal contri- 
butions come from the surface terms of the loop integrals they are suppressed in a renormalizable 
theory when A is kept large. Thus the gradient expansion anzatz can be justified for the evolution 
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kdkU{(P) 



IGvr' 



:U^'\<P). 



(42) 



The evolutions ( ^T] ) and P^ ) agree up to an overall constant in the scale parameter if 



dm'^ 



In 



'm^ + A^' 



m^ 



(43) 



can be considered as constant, i.e. m,^ and k independent. We obtain in this manner the 
usual justification of calling the Callan-Symanzik equation a renormalization group method, 
the equivalence of the scales in the internal and the external spaces in the ultraviolet scaling 
regime m'^ ^ U^"^^ up to non-universal contributions. The equivalence of the scales and the 
elimination of the non-universal contributions requires that the cutoff should be far above 
the mass, m'^ <^ A^. 

The non-vanishing anomalous dimension sets up the relation between the internal and 
external space scaling. In fact, when Z 7^ 1 the relation (^) becomes field dependent 
according to the first equation of (^). It is worthwhile comparing what (^) gives in the 
asymptotical regime m^ ^ f/*-^\ 



^m^ ^ni^ \ 



1 



32vr^Z^. 



In 



'Z„2(0)A2 + m^ 



m" 



z„A^)z''^A4>) 



[z^SM) 



(44) 



with the prediction of the Wegner-Houghton equation. A possible attempt to save the 
gradient expansion with sharp cutoff for the latter is the following: The contributions to 
the coefficient functions of the gradient, such as Zk{4>), come from taking the derivative of 
the loop integral, the trace in the second equation of (|^), with respect to the momentum of 
the infrared background field 0(x). There are two kind of contributions, one which comes 
form the derivative of the integrand, another from the external momentum dependence of 
the limit of the integration. It is easy to verify that the e-dependent non-local contributions 
come form the second types only [Q. Thus one may consider the approximation where 
these contributions are simply neglected, assuming a cancellation mechanism between the 
successive blocking steps. The result is, for k'^ ^ Ul (0), c.f. Appendix C, 



kdkZk{4>) 



k^ 



32vr2Z|(0) 



2Z,(</))zf(0) 



{Z^\^) 



(45) 



The formal similarity between the two different schemes, ( |4^ ) and (p5D, can be considered 
as a measure of the cancellation of the non-local terms evoked above. 
The beta-functions of the coupling constants Qn and Zm introduced as 



w) = eV' ^(<^)=eV 



n\ 



(46) 



are of the form 



I3n 



rn?drr,2 



Cdm 



m^9n 
Qn 



dyy^' 



Z(0)y + m2 + f/(2)(0) 



12 



7n = m^dm^Zn 



Qn „a2 
9</)" Jo 



Z(2)(0) 



+ 2Z^^\<P) 









2Z(0) 



Z«((/))y + f/(3) 



(Z(0)y + m2 + [/(2)(0))' 









(47) 



with Cd = i7d/2(27r)'^. The integration over y produces simple expressions for (3n and 7,1 
in terms of the couphng constants Qm and Zm- The simultaneous integration of this set of 
equations produces the solution of the evolution equation (P^). 

It is instructive to consider the solution in the independent mode approximation where 
the w? dependence is ignored in the integrals, C/((/)) = Ub{4>) ^^^ Zi.4>) = 1- We get 

f/e,,(0) = f/B(0) + I f dm' I / 



Us{<j^) + IJHp' + u^^\<P)] + o{m- 



which reproduces the usual one-loop effective potential for M ^ A. For the kinetic term, 
the integration of (B^) in the same approximation leads to 



Zeffi^) = 1 - 



1927r2 7m2 



dm 



K(0) 



m^ + U^^\^) 



uS\<P) 



1927r2 



uS\<P) 



0{M 



-2\ 



(48) 



for d = 4: which reproduces the one- loop solution found in [|1J]. The agreement between the 
independent mode approximation to the internal-space renormalization group equation and 
the one-loop solution is expected not only because the right hand side of (|29|) is 0{h) but 
because the one- loop contribution to the gamma function is universal, scheme independent. 
But this agreement does not hold beyond 0{h) as indicated by the imcompatibility of (^if ) 
and 



V. h DEPENDENCE 

It may happen that the quadratic suppression is not well suited to a problem. In the 
case 5*5 [0] possesses local symmetries which should be preserved then the choice (^) is more 
appropriate. The application of our procedure for a gauge model can for example be based 
on the choice 
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Sb[A] = -^J dxF;,F^^'^ + S,f[Al 
Ss[A] : 



1 



4^1 



dxF'^^F^"'^ 



(49) 



where Sgf contains the gauge fixing terms and on the apphcation of a gauge invariant 
regularization scheme. As mentioned after eq. (0) we need a regulator to start with in order 
to follow the dependence on the amplitude of the fiuctuations. One may use lattice, analytic 
(asymptotically free models) or Pauli-Villars (QED) regulator to render p^) well defined. 
The explicit gauge invariance of Sg [A\ which was achieved by suppressing the gauge covariant 
field strength instead of the gauge field itself makes obvious the independence of the resulting 
fiow for the gauge invariant part of the action from the choice of the gauge, Sgf. When a 
non-trivial saddle point appears in the blocking step then it may develop a discontinuous 



evolution. The choice (^Sj) makes the saddle point approximatively "renormalization group 
invariant" . 

We present the evolution equation for the 0^ model with quartic supression. 



S 



B\ 



s. 



-idf^cpx 



,2 I 92 ,2 I 93 ,3 I 94 ,4 



The similarity of this scheme with the loop expansion suggests the replacement 

1 



1 

h 

which yields the evolution equation 



1+A=1+ 



9 



(50) 



(51) 



1 



dgT[<P] = --e-^[^-l^, 



9^ 



S_ 



e^[^-l + 4^. 



9^ 



The integration of the evolution equation from g^n = to gfm '■ 
resummation of the loop expansion, i.e. the integration between h^ 

1 gives 



(52) 



oo coresponds to the 



and h 



fin 



The gradient expansion anstatz ( p3| ) with Z 
9 y^ -'p 



(53) 



+ 



5'4 

+— 
4! 



30 / G{p) - / G{p,)G{p2)G{-p,-p2) {u^'\<P)+9-\93 + 94^ 

Jp Jp\,P2 

3( IG{p)y + 6<f)^jG{p) 



P1,P2,P3 



G{p,)G{p2)G{p,)G{-p, -P2-P3) (U^'H^) + 



9 ^94 



- 3 / G(j9i)G(p2)G(p3)G'(-pi - P2)G{~pi -p2- Ps) 

x{u^^\<P) + g-\g, + g,^)y 

- 40 / G{pi)G{p2)G{-pi - P2) iu^'\<P) + g-\g3 + g4 
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where we used the fact that the Fourier transform of the IPI amphtude for n > 3 and Z = 1 



IS 



/ e*(^i-^^+-+P"-^'")r(")(xi,---,x„) = (27r)'^<5(pi + ---+p„)[/(")(0)- (54) 

The propagator in the presence of the homogeneous background field cf) is given by 

9, -'-^ 



G{p) 



p' + U^'\<P)+9-'(p' + 92 + 9,4>+'-l^' 



(55) 



Since the momentum dependence in the right hand side of (|53|) is exphcit and simple the one, 
two and three loop integrals can be carried out easily by means of the standard methods. The 
successive derivatives of the resulting expression with respect to yield the renormalization 
group coefficient functions. 

The use of internal space renormalization discribed in this section shows that this method 
can be generalized to any kind of action Sg and not only to a quadratic suppression term, 
as shown in the previous sections. 

VI. SUMMARY 

The strategy of the renormalization group is developed further in this paper. Instead 
of following the evolution of the coupling constants corresponding to the same physics our 
renormalization group flow sweeps through models with different dynamics. The parame- 
ter of the flow is the scale of the quantum or thermal fluctuations. The result is an exact 
functional differential equation for the effective action. As a special case the functional gen- 
eralization of the Callan-Symanzik equation is recovered. A different choice of the "blocking 
transformation" allows us to control h and the resulting flow amounts to the resummation 
of the loop expansion. 

Our scheme can be considered as a renormalization group method in the internal space. 
The similarity of the renormalization group flow in the external and the internal space is 
shown for the local potential, the zero momentum piece of the effective action. The difference 
between the two schemes is the source of the anomalous dimension and it appears in the 
momentum dependent parts of the effective action. 

The novel feature of the method is its manifest gauge invariance. This is achieved by the 
possibility of characterizing the modes in the internal space while the gauge transformations 
are carried out in the external space. 

The integration of the internal space evolution equation provides an algorithm to solve 
models in a manner similar to the traditional renormalization group method. The only 
truncation is done in the gradient expansion of the effective action, in the Taylor expansion 
of the IPI functions in the momentum. We believe that this procedure is an interesting 
alternative to the stochastic solution method of the lattice regulated models. The drawback 
is that it is rather cumbersome, thought possible in principle, to increase the precision in 
the momentum and going to higher orders in the gradient expansion. The advantage is that 
it can be cast in inflnite, continuous space-time equipped with Minkowski metric. 
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Note added in proof: After this work has been completed we learned that a method 



presented for gauge models in ref. []15| is similar to ours in the case of mass dependence 
(section IV). |T^ gives a loop expanded solution of the exact equation, whereas our solution is 
built in the framework of the derivative expansion. Finally, our approach can be generalized 
to any kind of suppression action Ss which is compatible with the symetries as shown in 
section V. 



APPENDIX A: LEGENDRE TRANSFORMATION 

We collect in this Appendix the relations between the derivatives of the generator func- 
tional W[i] and r[0] used in obtaining the evolution equations for V. 
We start with the definitions 



W\j] + T[<p\ + XSs[<P\=j 



Wi'\ 



and 



The first derivative of F gives the inversion of (|A2|) , 

^x = Jx — ^Js,x- 

The second derivative is related to the propagator W^^'f'^^ = G^^ 



X2 



r(2) 



Sjx^ 



\c(2) =n-^ -A^^^) 

^'-^s,xi,X2 xi,a::2 s,xi,X2' 



^X2 



The third derivative is obtained by differentiating 



f(3) 

^ X\,X2,X-i 



'yi,y2,yz 
The inverted form of this equation is 



/ G^^ G^^ G^^ w'^'^^ —A 9*^^) 

/ 2,'i,J/i"-^X2,J/2 X'3,!/3 ?;i,S/2,J/3 ^'~^S,Xx,X2,Xz- 

•JV^ .V9..v:i. 



Ty(3) 



i,X2,X3 / ^^i'yi^^2'y2^X3,y3 [^ yi,y2,ys ~^ ^'^s,yi,y2,y3 

J VI .VO.'U'i 



' 1/1,3/2,5/3 

The further derivation gives 



f(4) 

^ xi,X2,x:j,X4, 



yi,y2,y3,y4,zi,z2 



^-1 /-y-l /-y-1 /-y~l U/(4) 

^ xi,yi^ X2,y2^ X3,y3^ X4,y4'^'^ yi,y2,y3,y4 



+G"^ G~^ W^^'> G~^ W^^^ G'^ G'^ 

' ^Xi,yi^X2,y2 2/1 ,»/2,^l 21,^2 22, 2/3 ,»/4 2^3,2/3 3:-'4, 

+G-^ G~^ ly(^) G~^ ly(^) G~^ G'^ 

' ^3,y3 X2,y2 2/3 ,2/2, 21 '-'21,22 22,2/1,2/4 2^1,2/1 2,'4, 



2/4 
2/4 



+G~^ G~^ W^^^ G~^ ly(^) G~^ G~^ 

I ^ xi,yi^ X4„yi"' yi,y4,zi^ zi,Z2 22,2/3,2/2 2^3,2/3 2-2,2/2 

-AS'(^) 

S,Xl,X2,X3,X4' 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 
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Its inversion expresses the four point connected Green function in terms of the IPI amph- 
tudes, 



Xl,X2,X3,X4 



J/1 ,2/2,^3,3/4,^1 ,22 



l,J/2,S/3,J/4 



^xi,yi^X2,y2^xs,y3^X4,y4 (^J- j/i,y2,y3,J/4 ~'~ s,yi, 

l^-^ yi,y2,zi ^ ^'-'s,yi,y2,zij '^zi,Z2 y- Z2,yj,m ^ ^"^ s,Z2, 

-(1 n Tr^^) -I- A ^^^^ \ n Tr^^) -i- a 'i^^) \ n n 

^Xi,y3^X2,y2 \^ y3,y2,zi ' '^'~^ s,y3,y2,zi J ^ zx,Z2 y- Z2,y\,y4 ^ s,Z2,yi,y4j '~^xi,yi'^X4, 

'^xi,yi'^X4„y4, y- yi,y4,zi ^ ■^'^ s,yi,y4,zi J ^ zi,Z2 \y 22,S/3,J/2 ^ ^"^ s,Z2, 



(A8) 



yz,y4,) ^X3,y3^x4,y4 



VA 



,y3,y2j ^^3,J/3^^2,!/2 



APPENDIX B: EVOLUTION EQUATION IN THE INTERNAL SPACE 

We give here some details on the computation of (|35D . To get the evolution equation 
of the potential part of the gradient expansion (^), one has to take a homogeneous field 
= 00 in (HD- But to distinguish the kinetic contribution from the potential one, a non 
homogeneous field 4>{x) = 0o + vi^) is needed, as well. Let k be the momentum where the 
field 7] is non- vanishing. Then the effective action can be written as 



m = V,U,(^o) + \jmv{-q) {z,(^o)q' + f/f (0o)) + Oiv',k 



fBll 



where Vd is the spatial volume. Thus we need the second derivative of the effective action 
in (pop up to the second order in fj to identify the different contributions. The terms 
independent of fj give the equation for Ux and the ones proportional to k'^fj'^ the equation 
for Zx- The contributions proportional to r/^ but independent of k yield an equation for f/} 
which must be consistent with the equation for Ux- The result is 



r(2) 



Pl,P2 



Zx{(Po)pl + uPi<Po)]S{pi+P2) 



[B2) 



^(1)/ 



r(3). 



+ / ffiq) Z\:\<P,){pi + q' + qp,) + f/^ (0o) S{pi + P2 + q) 

Jq L J 



91, Q2 
:x3 7A\ 



'(2), 



Zx'i(Pom + H + qiq2 + 2gipi) + f/fl^o) Kpi +P2 + qi + ^2) 



(4)/ 



+ 0{fi\k 



Finally one computes the inverse of the operator \M.p^,p2 + ^p'i,p2 ^^"^ expands it in powers 
of fj and k. The trace over pi and p2 needs the computations of terms like 

Tr{(pigi)(p2g2)i^(pi,P2)<^(pi +P2 + qi + q2)} = -j^iq^ ^ ^^"^ J ^^^(^'"P) (^3) 



and they lead to 



The consistency with the equation for f/j; is satisfied. 
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APPENDIX C: EVOLUTION EQUATION IN THE EXTERNAL SPACE 

The evolution for the blocking in the momentum space is given by the Wegner- Houghton 
equation ([10|). In order to obtain its simplified version in the gradient expansion one needs 
the trace of the logarithm of r^^-'„2- This time the trace has to be computed for \pi\ and 
\p2\ in the shell [k — 6k, k]. This implies that jpi + g] has to be in the shell, as well, because 
^n,P2 contains 5{pi + P2 + q)- This constrain implies the appearance of \^ in the equation 
for kdkZk{(j)). The terms containing ^/q^ are non-local and spoil the the gradient expansion. 
kdkZk{(f)) is proportional to the second derivative of the trace in (|T0|) with respect the 
momentum q of the infrared background field, 0. The trace can be written as a momentum 
internal over the shell [k — 6k, k]. There are contributions from the dependence on q of the 
integrand and the limit of integration. It is worthwhile noting that the later contains all 
non-local terms. 

There are two ways to rid the non-local contributions when the model is solved by the 
loop expansion, i.e. by means of loop integrals for momenta < j9 < A. One is to use lattice 
regularization where the periodicity in the Brillouin zone cancel the q dependence of the 
domain of the integration. Another way to eliminate the non-local terms is to remove the 
cutoff. Since the non-local contributions represent surface terms they vanish as A — >■ cxo. 

One may furthermore speculate that some of the non-local terms cancel between the 
consecutive steps of the blocking k ^ k — Ak for a suitable choice of the cutoff function 
/(k) in the propagator G^^^^p) = f {p / k)G~^ {p) . Ignoring simply the non-local terms the 
identification of the coefficients of the different powers of the gradient in the two sides of 
(p!0|) leads to 

2{2'kY ""l Zfe(0)P + f/f (0) 




kd,u,{<p,) = -,,,^ In I i^:r:.„ :%;::' \ (ci) 



d 



kdkZ, 



'k^kWo) 



hVtdk 



2{2tiY 






l^Z,(0o)A:2 + f/r^(0o) (z,(0o)A:2 + f/i 



^0. 



d {z,{<j>o)k^ + t/f (0o)) ^ {Zk{h)k^ + ui'\<j>o)) 

(z,(0o)A:2 + f/f (0o)) ^ {z,{<Po)k' + f/f (0o)) 

When k'^ :§> Ul (0) this gives P3| ) in dimension d = A. 
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